Splitting between quadrupole modes of dilute quantum gas in a two dimensional 

anisotropic trap 
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We consider quadrupole excitations of quasi-two dimen- 
sional interacting quantum gas in an anisotropic harmonic 
oscillator potential at zero temperature. Using the time- 
dependent variational approach, we calculate a few low- 
lying collective excitation frequencies of a two dimensional 
anisotropic Bose gas. Within the energy weighted sum- 
rule approach, we derive a general dispersion relation of 
two quadrupole excitations of a two dimensional deformed 
trapped quantum gas. This dispersion relation is valid for 
both statistics. We show that the quadrupole excitation fre- 
quencies obtained from both methods are exactly the same. 
Using this general dispersion relation, we also calculate the 
quadrupole frequencies of a two dimensional unpolarized 
Fermi gas in an anisotropic trap. For both cases, we ob- 
tain analytic expressions for the quadrupole frequencies and 
the splitting between them for arbitrary value of trap defor- 
mation. This splitting decreases with increasing interaction 
strength for both statistics. For a two dimensional anisotropic 
Fermi gas, the two quadrupole frequencies and the splitting 
between them become independent of the particle number 
within the Thomas-Fermi approach. 

PACS numbers:03.75.Fi, 05.30Jp 



I. INTRODUCTION 

In a series of experiments, Bose-Einstein condensates 
(BEG) have been produced by cooling a vapor of alkali 
atoms to a temperature of a few nanokelvin ||l| . This sys- 
tem opens up interesting perspectives in the field of many 
body physics. There has been much progress in the the- 
oretical understanding of this system 1^. In particular, 
the low-energy collective excitation spectrum of a Bose 
condensed dilute gas in a trap has been discussed analyt- 
ically by Stringari using a hydrodynamic approxima- 
tion and also by sum-rule approach. A few low-lying ex- 
citations have also been calculated analytically by using 
time-dependent variational approach Q . The low-energy 
excitation spectrum obtained by using time-dependent 
gaussian variational ansatz exactly coincides with the hy- 
drodynamic results in the limit of large particle number 
N. Also a similar type of scaling ansatz has been used to 
describe the time evolution of the condensate in the large 
N limit 1^-0] . Experimentally the low- lying collective ex- 
citation frequencies of a condensate have been measured 



both at zero temperature g and at finite temperature 
1^. These observed values of the collective oscillation 
frequencies are in agreement with theoretical results at 
zero temperature. 

After the discovery of BEG in alkali atomic gas, the 
behaviour of trapped Fermi gas is also in focus. It is also 
possible to trap the Fermionic atoms at very low tem- 
perature, where the quantum effects can be observed. 
There has been experimental progress towards cooling 
a Fermi gas into the degenerate regime (T < Tp) [^ . 
Several authors have studied the thermodynamic proper- 
ties In , c ollective excitation frequencies in the normal 
phase |1^"^15[| as well as in the superfiuid phase [ p^JT?) of 
a three dimensional trapped Fermi gas. 

The reduction in dimension of a quantum system is 
the subject of extensive studies in trapped Bose systems 
]l8|-p2[ as well as trapped Fermi systems With 
present technology one can freeze the motion of the 
trapped particles in one direction to create a quasi two 
dimensional quantum gas by tuning the frequency in the 
z direction. Recently, the lower dimensional BEG has 
been realized p^ , p5[ . 

In a deformed quantum system, angular momentum is 
not a good quantum number and the angular momen- 
tum states mix with each other. Also in the presence of 
a small deformation of trap, the states with angular mo- 
mentum quantum number +1 and ~l split. The partition 
function of a deformed two-dimensional harmonic oscil- 
lator is exactly the same as the partition function of a ro- 
tating harmonic oscillator, where the rotation frequency 
n is related to the trap anisotropy, f2 = {ujy—ujx)/2. Here 
ujx and ujy are the oscillator frequencies of the deformed 
oscillator. From this analogy we find that the splitting 
occurs between the +1 and —I angular momentum states 
and it is proportional to the trap anisotropy. Similarly for 
a quantum gas in a deformed trap, the degenerate mul- 
tipole modes with angular momentum quantum number 
4-^ and —I also split. The dipole excitation frequencies of 
a quantum gas in deformed harmonic trap are the trap 
frequencies along x and y directions, lOx and ujy. The 
frequencies and the splitting of the dipole mode neither 
depends on the statistics of the trapped gas nor on the 
interaction between the particles. In circular symmet- 
ric trap two quadrupole frequencies are degenerate and 
the degeneracy is lifted up by elliptic trap deformation. 
In the present work we consider quadrupole modes of 
quantum gas in a quasi two dimensional anisotropic trap. 
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We calculate the frequencies of quadrupole like modes 
for both bosons and fermions in a deformed harmonic 
trap analytically. For both statistics, we study the effect 
of interaction on the splitting between the quadrupole 
modes for arbitrary deformation of the trap. First we 
consider the case of a two dimensional deformed trapped 
Bose gas at zero temperature. We analyze the nature 
of small oscillations of the confined gas within the time- 
dependent variational method. In case of deformed trap 
the quadrupole mode couples with the monopole mode 
whereas the scissors mode remains decoupled. This is 
due to fact that the Hamiltonian is invariant under re- 
flection. We calculate the frequencies of the mixed type 
of modes using variational technique. Next we construct 
the most general type of excitation operator which cou- 
ples the quadrupole and monopole mode. Then we con- 
sider only the minimum energy excitation in sum rule 
approach. The excitation energy obtained from the sum- 
rule approach exactly matches with that obtained from 
the time-dependent variational technique even for an ar- 
bitrary number of particles. In the case of a symmetric 
trap, these modes can be identified as a monopole mode 
and a quadrupole mode. There is a possibility of ex- 
perimental verification of our results since the quasi two 
dimensional Bose condensed has been realized in MIT 

il- 

This sum-rule method allows us to calculate the 
monopole mode and the quadrupole modes of a Fermi 
system in a deformed trap where the equation of motion 
technique does not hold. These modes can be observed in 
nano structures like quantum dots. When the two-body 
interaction is long ranged, the sum rule method can be 
generalised. 

There has been no systematic theoretical study on 
the collective excitations of a two dimensional deformed 
trapped quantum gas at zero temperature. The purpose 
of this paper is to give an analytic description of the 
quadrupole excitation frequencies of a two dimensional 
deformed trapped quantum gas at zero temperature and 
to calculate the splitting of the quadrupole modes for an 
arbitrary deformation of the trap. 

The paper is organised as follows. In Sec. II, we model 
the quasi-two dimensional trapped Bose system. Using 
the time-dependent variational method we calculate the 
monopole and quadrupole excitation frequencies of a two 
dimensional deformed trapped interacting Bose gas. In 
Sec. Ill, using the sum-rule approach we derive a general 
dispersion relation for the quadrupole excitation frequen- 
cies of a two dimensional trapped quantum system inter- 
acting through the two-body potential. This relation is 
valid for both Fermi and Bose statistics. We apply this 
general dispersion relation to calculate the same excita- 
tion frequencies of a trapped Bose system. We show that 
the quadrupole excitation frequencies obtained from both 
methods are exactly the same. In Sec. IV, we consider 
trapped unpolarized fermions and apply the dispersion 



relation obtained from sum-rule approach to calculate 
the frequencies of quadrupole modes. In Sec. IV, we 
present the summary and conclusions of our work. 



II. COLLECTIVE LOW-ENERGY EXCITATION 
FREQUENCIES OF A TWO DIMENSIONAL 
DEFORMED TRAPPED BOSE GAS 

In BEC experiments, the trap potential can be ap- 
proximated by an effective three dimensional harmonic 
oscillator potential, with tunable trap frequencies uJz in 
the axial (z) direction and lOx , i^y in the transverse (x-y) 
plane. The alkali-metal vapors used in experiment are 
very dilute and the interparticale interaction is well de- 
scribed by the short range pseudopotential and the inter- 
action strength is determined by s-wave scattering length 
a. Here we consider the case when the interparticle inter- 
action is strongly repulsive. The Gross -Pitaevskii(GP) 
p6| energy functional of the trapped boson of mass m is 
given by, 



Em = 



J 2m 



(1) 
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+ ^i^x^ +^yy +^zZW{r,z)\ 



where, g 



4-iran'' 



r is the position vector in x — y plane 



and ipi'f'j z) is the condensate wave function. 

It has been shown by Baganato et al. |^ that for 
an ideal two dimensional Bose gas under harmonic trap, 
a macroscopic occupation of the ground state can ex- 
ist at temperature T < = ^/N/({2)huj/kB. With 
present technology it is possible to freeze the motion of 
the trapped particles in one direction to create a quasi- 
two dimensional Bose gas. In the frozen direction the par- 
ticles execute zero point motion. To achieve this quasi- 
two dimensional system, the frequency in the frozen di- 
rection should be much larger than the frequency in the 
x-y plane and the mean interactions between the par- 
ticles. Alternatively, the trap frequencies are such that 
HlUz ^ fJ- > Tt-ij-iq and ksT <^ fitOz, where is the chemical 
potential of the two dimensional Bose gas. 

For a quasi two dimensional system we may assume 
that the wave function in the z-direction is seperable and 
is given by. 
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(2) 



where Qz — \/l\JmLOz is the oscillator length in the z- 
direction. Now we integrate out the z-component in the 
three dimensional GP energy functional, then we get the 
effective energy functional in two dimension: 
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E 



d\[—\vm\' 

2m 



(3) 



+ f / rfV5'(r-r')|V'(r')|2|V'(rOpi 



where §2 = 2\/2'KhLOzaza is the effective coupHng 
strength in two dimension, a is the s-wave scattering 
length in three dimension and N is the total number of 
particles in the condensate. The same effective coupling 
constant is obtained in Ref . [|l9| . The effective interaction 
in two dimension is given by, 



(4) 



The chemical potential of quasi two dimensional Bose 



condensed state is /i = hujoy \/2TrNa/az. Recently, the 
two dimensional Bose condensed state has been realized 
in MIT |2^. In this experimental set up, they have 
loaded iV 10^ number of sodium atoms in a trap with 
trap frequencies u/z/^n — 790 Hz, wq/Stt = ^ujxiOy ~ 
20Hz and a = 2.75 nm. One can easily calculate the 
chemical potential fi ^ 0.19hLJo^/N which satisfies the 
above mentioned inequality condition to be a quasi two 
dimensional Bose system. 

In two dimensions, the equation of motion of the 
condensate wave function is described by the Gross- 
Pitaevskii equation. 



5V'(r) 



ih 



dt 



= h^V^ + V{r) + 52|V'(^1P]^(^"*), (5) 



where ^(r) — ^■m{ujlx^ + (jJyU^) is the deformed trap 
potential in two dimension. The normalization condition 
for ij} \s J d^rji/jp = TV. N is the number of particles in 
the condensate . One can write down the Lagrangian 
density corresponding to this system as follows: 



C 



'2m 



dt' 



(6) 
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where * denotes complex conjugation. One can get the 
non-linear Schroedinger equation (|^) by minimizing the 
action related to the above Lagrangian density (H). In 
order to obtain the evolution of the condensate we assume 
the most general Gaussian wave function. 



(7) 



where C(t) is the normalization constant. X and Y are 
the dimensionless variables, X = Y = where 

' ao ' ao 

flo \J'^^^ oscillator length and uoq — ^JUTJU^ is 

the mean frequency. Further, a = ai + ia2, /3 = Pi+ i(32 
and 7 = 7i + «72 are the time dependent dimensionless 



complex variational parameters. The ai and /3i are in- 
verse square of the condensate widths in x and y direction 
respectively. The square of the normalization constant is 
|C(t)|2 = i!^, where D = ai(3i - 7^2. The Gaussian 
ansatz (eq. ^) for the order parameter can also be gen- 
eralized to three dimensional anisotropic trapped Bose 
system to study various scissors modes. 

The gaussian variational ansatz becomes an exact 
ground state in the non interacting limit and in the 
presence of repulsive interaction it gives rise to spread- 
ing of the condensate wave function. To describe the 
quadrupoles and monopole oscillation, we consider the 
most general time dependent quadratic exponent of the 
variational ansatz. 

We obtain the effective Lagrangian L by substituting 
Eq. into Eq. (||) and integrating the Lagrangian 
density over the space co-ordinates. 



■\-{Pid2 + q;i/32 - 27172) 
-^12). 

Oil 



(8) 



+ (ai +(3i)D + [al + -fl)(3i + {Pi + 72)^1 
- 2(a2 + /32)7i72 + A/3i + ^ + PD''/^ 

where A is the asymmetric ratio, A = ojx/ojy and P = 
^2N-S-. 

The variational energy of the static condensate at equi- 
librium is given in terms of the equilibrium values of the 
inverse square width of the condensate along x and y 
directions, 



E 1 



A 



1 



(9) 



One can get the equilibrium value of the variational pa- 
rameters, aio and Piq by minimizing the energy with 
respect to the variational parameters, 



P 



T ^ irAovT^ioOi 



From the above two relations, we obtain 



J] 



Prf 



P-q 
2A2 



1 

A2 



0, 



(10) 



(11) 



(12) 



where is the ratio of the condensate widths in the x 
and y direction, 77 = -y/f^- From Eq. ( |l2|) one can 
say how 77 changes with the number of atoms N and the 
coupling constant 52- The variation of rj with the the 
dimensionless effective interaction strength P is shown in 
Fig. 1. The ratio between the widths of the condensate 
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r] varies from 1/VX to 1/A, as the interaction strength 
increases from zero to large value (Thomas- Fermi limit). 
In the Thomas-Fermi limit, the equilibrium values of the 
parameters ai and Pi are, 



(13) 



"10 — -^y pi Pio ^ ^ Y p 



In this limit, the energy per particle is ^ = hcjQy —. In 

the non-interacting limit, o^q = A and /S^q = 1/A. The 
energy per particle is ^ = -^(-v/A-I- a/1/A). 

We are interested in the low-energy excitations of a 
Bose system. The low-energy excitations of the con- 
densate correspond to the small oscillations of the cloud 
around the equilibrium configuration. Therefore we ex- 
pand the time dependent variational parameters around 
the equilibrium points in the following way, ai = aio + 
Sai, Pi = Pw + SPi and a2 = Sa2,P2 = SP2 71 = ^71, 
710 = and 72 = (572. 

Using the Euler-Lagrange equation, the time evolution 
of the inverse square of the width around the equilibrium 
points are given by, 



Sdi + A- 



(2 + P77) 



dai 



PXri 



(2 + Pr^) 



SPi = 0, 



(14) 



SPi + 



Xi2r] + P) 



(877 + 3F),^ ^ 



XiP + 2r]) 



^7r + [(^ + (A + >7i=0. 



(16) 



From the above equations we can see that the modes 
corresponding to the fluctuations of the average of x"^ 
and are coupled, but the mode associated with the 
fluctuation of the average value of xy is decoupled. The 
Eqs. ( |l^ ) and ( 15) are coupled equations of the modes 
ai and Pi, where the mode xy is decoupled. This is 
due to the fact that the Hamiltonian is invariant under 
reflection, x —x and y ^ y or x ^ x and y —y, 
and the modes which are odd or even under this operation 
separate out. 

Now we look for time dependent solutions of e*'^* type, 
we obtain from Eqs. ( [T^ ) and (p^. 



X{8 + 3Ptj) 1 (877 + 3P) 
2 {2 + Pf]) ^ 2A (277 -HP) 



(17) 



^ A(8-f-3P?7) 

2(2-fP?7) 2A(2r7 + P)^ 



+ (7 



Pf]' 



For an isotropic trap, lu-^ — 2ujq and uj"^ — [wq(8 + 
2P)]/(2 + P). For large N limit, = \/2ujq. So 
loj^ and lo- may be identified as the monopole mode 
frequency and quadrupole mode frequency respectively. 



The monopole mode is coupled with the quadrupole 
mode in an anisotropic trap. However, the monopole 
mode frequency in an isotropic trap is independent of 
the interaction strength of the two-body potential and 
the number of particles in the condensate state. This is 
due to the underlying SO{2, 1) symmetry in the Hamil- 
tonian |§, pi. 

From Eq.(16) we obtain. 



Lot 



AXrf 



(2 + Pr^) 



+ (A + i) 



(18) 



In an isotropic trap, lJs becomes lo-. In the non- 
interacting limit, ujs = LUx + ujy. In the Thomas-Fermi 

limit, Eq. ( p^ ) reduces to lUs ~ \J^x + ■ 
isotropic trap this mode corresponds to the quadrupole 
excitation. This excitation is also known as scissors mode 
p8| , and this oscillation has been observed experimen- 
tally M. 



III. SUM-RULES AND COLLECTIVE 
EXCITATIONS 

In this section, we study the quadrupole excitations 
of a two dimensional deformed trapped quantum gas at 
zero temperature within the sum rule approach. In the 
collisionless regime the collective excitation frequencies of 
a confined gas are well described by the sum rule method. 
The collective excitation of any system is usually probed 
by applying external fields. Given an excitation operator 
P, many useful quantities of the excited system can be 
calculated from the so-called strength function jsol , 



^±(£^) = E I < "l^±|o > ^»)' 



(19) 



where En and \n > are the excitation energy and the ex- 
cited state respectively, and P+ = F, F- = F^ . Various 
energy weighted sum rules are derived from the moments 
of the strength distribution function. 



i I E''iS+iE)±S^iE))dE. 



(20) 



It is easy to see that, for a given fc, the moments may be 
expressed in terms of the commutators of the excitation 
operator F with the many body Hamiltonian H. We give 
below some of the useful energy weighted sum rules, 



^0 -^<oi[pnp]io>. 



-<0|[pt,[i7,P]]|0>, 



<0|[J^J]|0>, 



(21) 



(22) 



(23) 
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4 ^^<mF\H],[H, [H,F]]]\0>,J 



[H,F], (24) 



where [, ] denotes the commutator between corresponding 
operators. Near the collective excitation frequency the 
strength distribution becomes sharply peaked, and the 
collective excitation energy is described by the moments 
of the strength distribution, 




(25) 



Following Ref. |31|, one can derive the above form of 
collective excitation energy by using the variational prin- 
ciple. Given the many body ground state it is possible to 
find out the collective excitation energy and the excited 
state, if one is able to find an operator , which satisfies 
the following equation of motion : 



(26) 



The excitation energy is then given by the following ex- 
pression. 



huJr ~ 



<0|[0,[g,Qt]]|o> 
<0|0,0t|o> ■ 



(27) 



We may now take the variational ansatz for as, 
= F + bJ with the variational parameter b. By min- 
imizing the energy with respect to the variational pa- 
rameter, we obtain the collective excitation energy as 

Ec — \J /m^ ^ which is same as eqn.(p5|). 

Similarly, we construct the most general excitation op- 
erator F = -\- by^ when monopole and quadrupole 
modes are coupled. & is a variational parameter. In sym- 
metric trap potential, if & = 1, F is monopole mode and 
if 6 = —1, F is the quadrupole mode. In the same way we 
can calculate the lowest energy excitation in this particu- 
lar sector of excitations variationally. The lowest energy 
mode turns out to be the quadrupole mode. 

Calculating the moments mi and ms by taking the 
excitation operator with the Hamiltonian, we obtain. 



^coll — 



Ah' 



E^ (1 + Ab^ + Bb) 



m < > {1 + Ch^) 



where A = 



2E^ 



,c 



'J < Eint >, Ey 



<y > 

--<Ty > 



and E.r 



(28) 



=< > 



<Vy> 



< 



Eint >■ Here <> denotes the expectation value of the 
corresponding operators in the ground state and Tx, Vx 
and Ty, Vy represents the kinetic energy and potential 
energy along x and y coordinates respectively. The in- 
teraction energy is given by Eint — 52/2 / \^p\^d'^r. Now 
we minimize this collective energy with respect to the 
variational parameter p. The value of /3 for which the 
collective energy is minimum, is given by 



60 



-2{C -A)± ^A{C - A)"^ +4B'^C 
2BC ■ 



(29) 



It can be easily shown that for an isotropic trap, bo = ±1. 
So we have identified that the monopole mode can be ex- 
cited by the operator F — x'^ + b^y^ and quadrupole 
mode can be generated hy F = x^ — boy^. Inserting bo 
into Eq. (p8|), we obtain the following collective oscilla- 
tion frequencies: 



E,i 



Ex 



-'± — LV 9 

m < y 



> 



< x^ > 



(30) 



± 



E„ 



Ex 



E^ 



■< r > 



<x^ > 



A<x^ >< y^ > 



Using the variational wave function of the ground state 
in deformed trap, one can easily get the excitation fre- 
quencies. The lowest energy excitation frequency in this 
sector is 



X{8 + 3Pr)) 1 {8ri + 3P) 
2 (2-|-P?7) ^ 2A (2?7 + P) 



(31) 



r A(8 + 3Pr;) _ (8r? + 3P) ^ ^ PXy' 



Pf]' 



It can be identified as quadrupole mode since in an 
isotropic trap, its excitation frequency exactly matches 
with the quadrupole mode frequency. The above expres- 
sion for the excitation frequency Eq. ( ^l| ) is exactly same 
as the mode frequency to- in Eq. (p^). 

The higher energy excitation exactly matches within 
the monopole mode, although it is not the local minimum 
of the energy, Eq. ( p8| ) . The dispersion relation for this 
monopole mode frequency is: 



X {8 + 3Pr]) J_(877 +_3P) 
2 (2 + Pr/) ^ 2X (2?7 + P) 



(32) 



.A(8 + 3Pr;) (877 + 3P) 



PA772 



2(2-|-P77) 2A(277 + P) 



Pri' 



For another quadrupole mode, the excitation operator 
is P = xy. Using the commutation relation, we obtain, 



mi 



2m 



< (a:2 + y2) >^ 



where, 



^3{T) = — < (pi+p;) >, 



ms{V) = —{ui+u;'y)<{x'+y')> 



(33) 



(34) 



(35) 



(36) 
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52^ 

2m2 



(frp{r){y 
dp 



dy 



y 



d_ 

dx 
dp, 21 



+ ^-^fpir) (37) 



dx' 



Using the variational wave function, wc can get all 
these moments. In this case m3(ee) exactly vanishes. 
So the frequency for the quadrupole mode is, 



LU^ (2 + Pr,) 



(A- 



1. 



(38) 



This expression for the quadrupole frequency is also the 
same as Eq.(|l|). So w_ in Eq. @ and w,. in Eq.(||) 
shows the splitting occurs between two quadrupole modes 
in a two dimensional deformed trapped Bose gas and the 
dependence of the splitting on the interaction strength 
and trap anisotropy can be analyzed from the analytical 
expressions. For an isotropic trap, the two quadrupole 
modes are degenerate. The variation of the splitting be- 
tween two quadrupole modes Af, = — oj-, of a trapped 
interacting Bose gas, with the dimensionless interaction 
parameter P is shown in Fig. 2. 

We have checked that the sum rule method gives cor- 
rect results for the excitation frequencies of the two 
quadrupole modes for a system of interacting bosons in 
an anisotropic trap. Now we apply this method to calcu- 
late the excitation energies of the quadrupole modes of a 
system of interacting Fermions in a deformed trap. 



The Thomas-Fermi energy functional of this trapped 
interacting Fermi system is given by, 

E[p{r)] = I d'ri^p'ir) + V,,,p{r) + ^p\r% (41) 

where 172 — 52/2. Here we assume the density of two spin 
components are same, p| — Pi- The interaction energy 
density g2PiPi, can be written as ^p^, where p is the 
total density. By minimizing the energy functional with 
respect to density, we obtain. 



p(r) 



],r<RF, 



(42) 



where Rp = {ANKoY^^ao is the radius of the atomic 
gas which is determined from the condition J <Prp{r) = 
N . Kq = 1 + is a dimensionless constant. At very 
low temperatures, collisions are suppressed due to Fermi 
statistics and system is in the coUisionless regime. We 
study the collective excitation frequencies in this regime 
by sum rule approach. 

In Sec. Ill we have derived the expressions for 
quadrupole excitation frequencies within sum-rule ap- 
proach. We can use expression Eq.(pC|), to calculate one 
of the quadrupole mode frequencies for Fermi gas also. 

We evaluate all the expectation values of the corre- 
sponding operators by using the Thomas-Fermi density 

(1), 



IV. TWO DIMENSIONAL TRAPPED 
ANISOTROPIC FERMI GAS AT ZERO 
TEMPERATURE 



< y' > 2 2ao a 



Using Eqs. (|30|) and (43), we obtain, 



In this section we discuss the collective oscillation of 
a two dimensional deformed trapped unpolarized Fermi 
gas at zero temperature within the sum-rule approach. 
Using this approach, the collective excitations have been 
studied in other finite Fermionic systems like atomic nu- 
clei, metal clusters and quantum dots p3[ |. 

We consider a two dimensional deformed trapped un- 
polarized Fermionic atoms at very low temperature. The 
two-body interaction of the dilute gas can be described by 
the short range pseudopotential V{f— f) = g26^{f— f'), 
where g2 is the coupling constant and its form is given in 
Sec. II. The Hamiltonian of the trapped Fermionic atoms 
is given by. 



H = 



PL 



g2y^J^in-rj), 



where the confining potential is 



1 



Vext = -mwo(Ax + Y) 



(39) 



(40) 



(3X, 



(44) 



± 



For an isotropic trap, the monopole mode frequency be- 
comes — 2luo. 

There is another quadrupole mode for which the ex- 
citation operator is F = xy. Using the density for the 
trapped interacting Fermi gas at T = 0, we get the fol- 
lowing moments: 



■nil 



n^R% 



(A 



12A'om2ag'Xo 



In this case also, m3(ee) exactly vanishes, 
quadrupole oscillation frequency is given by. 



(45) 

(46) 
The 
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For an isotropic trap, ui- in Eq. (Q) and Eq. ( ^7| ) be- 
comes 



.., = 72^.0^^^^. (48) 

In Eq.(|^ and Eq.(|^), lj_ exhibits the spUtting of 
the quadrupole modes of a two dimensional deformed 
trapped Fermi gas. The scissors mode is also discussed 
in Rcf. for superfluid Fermi gas. 

The monopole mode frequency of an isotropic trapped 
interacting Fermi system is 2u]q which is independent of 
the interaction strength of the two-body potential and 
the number of particles. This is because of the presence 
of S0(2,l) symmetry in the Hamiltonian (Eq.(p9[)) as 
discussed in the previous section. 

The splitting between two quadrupole modes, Aj^ = 
ujs — w_ , of a deformed trapped interacting Fermi gas is 
shown in Fig. 3. The frequencies of these two modes and 
the splitting between them are independent of the parti- 
cle number for two dimensional Fermions within Thomas- 
Fermi approximation. This splitting decreases almost lin- 
early with increasing interaction strength. 

V. SUMMARY AND CONCLUSIONS 

In this paper, we have mainly considered two non- 
degenerate quadrupole modes of a quantum gas in an 
anisotropic harmonic oscillator potential. We investi- 
gated the effect of interaction on the splitting between 
these quadrupole modes for arbitrary trap deformation. 
We have calculated a few low-lying collective excitation 
frequencies of a two dimensional trapped Bose gas in 
an anisotropic trap, by using time dependent variational 
method. We found that one quadrupole mode is coupled 
with the monopole mode in presence of trap deformation. 
Another quadrupole mode associated with the fluctua- 
tion of the average value of xy (which is also known as 
scissors mode), is decoupled. 

Using the energy weighted sum-rule approach we 
derived the general dispersion relation of the two 
quadrupole excitations. Using the same variational wave 
function for Bosons, we checked that the collective fre- 
quencies obtained from the sum-rule approach are ex- 
actly the same as those obtained from the variational 
method. The main advantage of the sum-rule method 
is that it can be applied to both trapped Bosons and 
Fermions to calculate the excitation frequencies in the 
coUisionless regime. This method can be applied for 
any number of confined particles and also it can be 
generalised for long range interactions. This energy 
weighted sum-rule method can be extended for Coulomb 



interaction to study the quadrupole excitations in a de- 
formed electronic nanostructure like an elliptic quantum 
dot. The splitting between the quadrupole modes ob- 
tained from this method is non perturbative in the trap 
anisotropy parameter. 

We considered a system of two dimensional spin un- 
polarised interacting Fermions in an anisotropic har- 
monic oscillator potential within Thomas-Fermi approx- 
imation. Applying the sum-rule technique to this de- 
formed Fermi gas, we obtain two quadrupole excitation 
frequencies and the splitting between them analytically. 
For both statistics, the amount of splitting between the 
quadrupole modes decreases with increasing interaction 
strength. For a two dimensional Fermi system the fre- 
quencies and the splitting are independent of the parti- 
cle number. For an isotropic trap, the monopole mode 
frequency of a Bose gas as well as Fermi gas is the uni- 
versal frequency 2lO{). This monopole mode frequency is 
independent of the strength of the two-body interaction 
potential and the number of particles. This is due to 
the underlying 50(2, 1) symmetry in the Hamiltonian. 
Strictly speaking, our all the results are valid when the 
conditions Tiojz ^ > Uloq and ksT <^ huz are satisfied. 

Recent experimental progress in MIT |Q on quasi two 
dimensional Bose condensed shows the possibilities of 
verification of our results. Above mentioned quadrupole 
modes are excited in the two dimensional plane and for 
simplicity we consider only the two dimensional trapped 
gas. This method and the most general Gaussian anstaz 
for the order parameter can also be extended to three 
dimensional anisotropic systems to study the various 
quadrupole modes. The splitting in these two quadrupole 
modes may be used to find trap anisotropy. It will be an 
interesting to study the splitting between the quadrupole 
modes of an anisotropic quantum system in presence of 
terms having definite chirality, like magnetic field or ro- 
tation. 
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FIG. 1. The variation of the ratio of the widths of the 
condensate, ry as a function of the dimensionless effective in- 
teraction strength P, for the fixed ratio of trap frequencies 
A = 0.7. 

FIG. 2. The difference between the two quadrupole modes 
of an interacting Bose gas, Ab/uJo as a function of the dimen- 
sionless effective interaction strength P for fixed ratio of trap 
frequencies A = 0.7. 

FIG. 3. The difference between the two quadrupole modes 
of an interacting unpolarized Fermi gas Af/uo, as a function 
of the dimensionless parameter Ko for fixed ratio of trap fre- 
quencies A = 0.7. 
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